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I.  Introduction 


This  project  is  concerned  with  the  fabrication  and  investigation  of  superconducting  flux  qubits, 
particularly  with  regard  to  the  effects  of  dissipation  on  their  relaxation  and  decoherence,  and  to  the 
entanglement  of  two  qubits.  This  final  report  summarizes  the  progress  that  we  have  made  towards 
achieving  these  goals  over  the  last  three  years.  Section  II  outlines  the  infrastructure  of  our  experiment, 
including  the  dilution  refrigerator,  the  electrical  filtering  of  the  various  leads  associated  with  the 
measurements,  a  novel,  dual-current  source  to  supply  two  independent  flux  biases  to  the  qubits  and  readout 
SQUID,  the  chip  layout,  and  the  fabrication  process.  Section  III  describes  our  results  on  achieving 
quantum  coherence  in  a  flux  qubit,  including  spectroscopy,  Rabi  oscillations,  Ramsey  fringes, 
spectroscopic  linewidths,  and  spin  echoes.  The  different  kinds  of  decoherence  times  deduced  from  the  last 
three  experiments  are  shown  to  be  self-consistent.  Section  IV  describes  a  study  of  the  effects  of 
nonequilibrium  quasiparticles  generated  in  the  readout  SQUID  when  it  switches  to  the  voltage  state.  These 
excess  quasiparticles  are  shown  to  persist  for  a  remarkably  long  time,  about  1  ms,  and  thus  set  an  upper 
limit  on  the  repetition  rate  at  which  data  can  be  acquired.  Section  V  describes  the  theory  of  a  novel  device 
for  entangling  two  flux  qubits  by  means  of  a  single  SQUID  that  serves  also  as  the  readout  device.  It  is 
shown  that  this  scheme  is  in  principle  capable  of  the  CNOT  (Controlled  NOT)  operation.  Section  VI 
contains  our  concluding  remarks. 

II.  Experimental  Configuration 

In  order  to  observe  quantum  coherence  in  a  superconducting  circuit,  it  is  necessary  not  only  to 
cool  it  to  low  temperatures— say,  a  few  tens  of  millikelvin — but  also  to  isolate  it  from  thermal  and 
environmental  noise  generated  at  room  temperature  and  at  lower  temperatures  within  the  refrigerator  itself. 
Thus,  one  needs  to  design  a  series  of  filters  and  attenuators  that  provide  large  levels  of  attenuation  at  room 
temperature,  4.2  K  and  certain  temperatures  between  4.2  K  and  the  temperature  of  the  experiment. 
Furthermore,  it  is  of  great  importance  that  the  qubit  and  its  readout  SQUID  (Superconducting  Quantum 
Interference  Device)  are  coupled  to  impedances  that  are  high  enough  to  ensure  that  any  Johnson  noise 
currents  injected  into  the  quantum  circuit  are  unimportant  yet  devoid  of  parasitic  resonances  at  frequencies 
relevant  to  the  operation  of  the  qubit,  typically  a  few  gigahertz. 
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Fig.  1  Experimental  configuration. 
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Figure  1  shows  schematically  the  configuration  of  the  measurement  system.  The  dilution 
refrigerator  is  surrounded  by  a  solid  copper  shielded  room  with  an  attenuation  of  100  dB  at  1  GHz.  The 
refrigerator,  an  Oxford  Kelvinox  300,  is  capable  of  attaining  temperatures  below  20  mK;  however,  the 
thermal  load  of  the  various  wires  and  cables  limits  the  base  temperature  to  35^40  mK.  The  chip  on  which 
the  superconducting  circuits  are  fabricated  is  mounted  in  the  copper  box.  The  cavity  surrounding  the  chip 
is  coated  with  lead  and  closed  with  a  lead-coated  lid  (not  shown)  so  that  the  circuits  are  enclosed  in  a 
superconducting  cavity.  The  copper  box  is  thermally  anchored  to  the  mixing  chamber.  The  pulse  and 
sense  lines  for  the  readout  SQUID  enter  from  the  right-hand  cavity,  and  are  heavily  damped  by  3  kO 
resistors.  The  two  flux  bias  lines,  one  for  the  qubit  and  the  other  for  the  readout  SQUID,  enter  from  the  left 
and  are  heavily  damped  by  copper  powder  filters  [  1]  in  the  cavity  (not  shown).  Finally,  microwaves  can  be 
coupled  to  the  qubit  via  a  coaxial  line  that  enters  from  below  and  is  coupled  to  the  qubit  via  a  1-mm- 
diameter  loop  of  niobium  wire. 

Figure  1  shows  the  various  levels  of  filters  and  attenuators.  To  determine  the  critical  current  of  the 
readout  SQUID,  a  current  pulse  is  supplied  by  a  computer-controlled  generator  outside  the  shielded  room. 
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This  line  has  a  20-dB  attenuator  at  room  temperature,  a  20-dB  attenuator  at  1 .5  K,  a  10-dB  attenuator  at  0.7 
K  and  a  6-dB  attenuator  at  40  mK.  The  pulse  is  largely  reflected  by  the  3-kQ  terminating  resistor,  but  the 
reflected  pulse  is  absorbed  by  the  matched  50 -ft  attenuator  at  40  mK  so  that  there  is  no  resonance.  The 
sense  line  is  coupled  to  the  readout  SQUID  via  a  3-kft  resistor,  and  the  signal  passes  through  a  copper 
powder  filter  at  40  mK  and  a  Tt-filter  at  1 .5  K  before  being  coupled  to  the  room-temperature,  battery- 
operated  preamplifier.  If  the  SQUID  switches  in  response  to  the  pulse,  the  comparator  transmits  a  signal  to 
the  computer  outside  the  shielded  room  via  a  fiber-optic  link.  The  flux  bias  currents  are  very  heavily 
filtered:  there  are  Tt-filters  at  300  K  and  4.2  K,  a  copper  powder  filter  at  4.2  K,  and  a  7i-filter  and  two 
copper  powder  filters  at  40  mK.  This  series  of  filters  roll  off  sharply  at  frequencies  above  25  Hz.  Finally, 
the  microwave  generator,  which  is  outside  the  shielded  room,  is  inductively  coupled  to  the  qubit  via  a  line 
with  the  following  attenuators:  20  dB  at  room  temperature,  10  dB  at  4.2  K,  3  dB  at  1.5  K,  3  dB  at  100  mK 
and  3  dB  at  40  mK.  In  addition,  there  is  a  high  pass  filter  at  10  MHz.  The  mutual  inductance  of  the  line  to 
the  qubit  is  extremely  small,  approximately  10  fH. 

The  low-temperature  stage  of  the  refrigerator  is  surrounded  by  a  cryogenic  mu-metal  shield  that 
reduces  the  earth’s  magnetic  field  to  below  1  pT. 


to  potentiometers  of  channel  2 

Fig.  2  Schematic  of  current  source. 


Flux  bias  currents  are  produced  by  digitally  programmable  potentiometers  controlled  by  computer 
over  a  fiber-optic  link  [2].  This  circuit,  shown  schematically  in  Fig.  2,  consists  of  two  current  channels  each 
with  a  maximum  output  of  ±1  mA  and  an  accuracy  of  about  ±1  nA.  The  current  source  was  designed  not  to 
inject  digital  noise  into  the  qubit,  and  measured  noise  currents  are  sufficiently  low  to  contribute  negligibly 
to  observed  decoherence.  The  low  noise  of  the  current  source  combined  with  the  shielding  from  the 
superconducting  cavity  provides  high  flux  stability  in  the  system,  allowing  for  experimental  runs  of  up  to 
48  hours. 


3 


Fig.  3  Chip  layout.  Blue  represents  A1  traces,  gold  AuCu  traces.  Pads  near  upper  edge  of  chip  provide 
two  independent  flux  lines;  A1  wirebonded  jumpers  couple  left  and  right  halves  of  each  of  the  two  flux 
lines.  Pads  near  lower  edge  of  chip  are  used  to  supply  current  pulses  to  the  read-out  SQUID  and  to  sense 
any  resulting  voltage.  Chip  is  2  mm  on  a  side,  (b)  Photograph  of  completed  device.  Segments  of  flux  lines 
are  visible  to  left  and  right  of  SQUID,  which  surrounds  the  two  qubits. 


Figure  3  shows  the  layout  of  our  device.  The  read-out  SQUID  has  a  calculated  inductance  Ls  = 

358  pH,  and  each  of  the  two  qubits  it  encloses  has  a  calculated  inductance  L0  =  143  pH.  The  calculated 
mutual  inductance  between  each  qubit  and  the  SQUID  is  6 1  pH.  One  pair  of  series-connected  flux  bias 
lines  is  arranged  near  the  top  of  the  SQU  ID  and  a  second  near  the  bottom;  thus,  flux  in  qubit  1  (2)  is 
supplied  predominantly  by  the  upper  (lower)  flux  lines.  The  mutual  inductance  between  each  qubit  and  its 
associated  flux  lines  was  designed  to  be  4-5  pH,  enabling  us  to  apply  ~2<D0  with  a  current  of  1  mA.  This 
criterion  dictated  the  relatively  large  qubit  inductances  compared  with  those  in  previous  experiments  [3,4]. 
The  wiring  to  the  SQUID  that  provides  current  pulses  and  senses  any  resulting  voltage  was  laid  out 
symmetrically  to  minimize  spurious  inductive  coupling  to  the  SQUID  and  qubits. 

We  fabricated  the  device  on  an  oxidized  Si  substrate  using  electron-beam  lithography  and  double¬ 
angle  evaporation  to  form  the  Al-AlOx-Al  tunnel  junctions.  The  Al  lines  for  the  qubit  and  SQUID  loops  are 
1  pm  wide,  and  for  the  flux  bias  1 0  pm  wide.  Each  SQUID  junction  was  1 75  x  200  nnr  with  a  capacitance 
of  6.5  fF,  estimated  with  a  specific  capacitance  of  100  fF/pnr  and  a  stray  capacitance  of  3  fF,  obtained  by 
separate  experiments  on  SQUIDs  with  a  range  of  areas.  The  junction  critical  current  was  0.44  pA.  For 
qubits  1  and  2,  the  larger  junctions  had  areas  of  250  x  250  nm2  and  1  SO  x  200  nnr,  critical  currents  I0 
(scaled  from  SQUID  junctions)  of  0.79  pA  and  0.46  pA,  and  a-vaiues  of  0.49  and  0.68,  respectively.  The 
different  junction  parameters  were  chosen  to  increase  the  probability  of  obtaining  one  viable  qubit;  in  fact, 
qubit  2  displayed  good  characteristics.  A  AuCu  film  deposited  and  patterned  prior  to  the  Al  deposition 
provided  quasiparticle  traps  near  the  junctions,  shunts  on  each  of  the  four  flux  lines,  and  series  resistors  on 
each  end  of  the  pulse  line  and  sense  lines. 
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III.  Quantum  Coherence  in  a  Flux  Qubit 
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Fig.  4 ,  Spectroscopy  of  flux  qubit. 


We  have  observed  quantum  coherence  in  qubit  2  [5].  Figure  4  shows  the  excitation  frequency 
versus  the  qubit  flux  obtained  in  a  computer-controlled  experiment  in  which  the  microwave  frequency  and 
qubit  flux  are  varied  with  a  constant  flux  in  the  readout  SQUID.  Spectroscopic  peaks  are  on  the  right  and 
dips  on  the  left;  there  is  a  well-defined  splitting  of  A/h  =  3.99  GHz  at  the  degeneracy  point.  The  expected 
cavity  resonance  is  clearly  visible  at  about  6  GHz.  From  this  spectrum,  we  deduce  a  flux-to  frequency 
conversion  of  896  MHz/mel>0l  The  outer  two  branches  correspond  to  two-photon  transitions,  fitted  with  the 
same  parameters  as  the  single-photon  spectrum. 
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Fig.  5.  (a)  Rabi  oscillations  of  flux  qubit  for  three  different  microwave  powers,  (b)  Rabi 

frequency  vs,  microwave  drive  amplitude. 


Figure  5(a)  shows  a  series  of  Rabi  oscillations  obtained  by  applying  a  pulse  of  microwaves  to  the 
qubit  just  before  measuring  its  flux  state.  The  switching  probability,  scaled  to  the  fidelity  of  the  readout 
SQUID,  is  plotted  versus  the  duration  of  the  microwave  pulse.  Fits  of  exponential  decaying  sine  waves  to 
the  data  reveal  relatively  large  contrasts  -  as  high  as  63%  -  and  Rabi  decay  times  as  high  as  78  ns.  Figure 
5(b)  shows  that  the  measure  Rabi  frequency  scales  linearly  with  the  microwave  amplitude,  as  predicted. 
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Fig.  6 ’  (a)  Ramsey  fringe  at  150  MHz  detuning,  (b)  Ramsey  fringe  measurement  containing  80 

excitation  frequencies. 


Figure  6  illustrates  Ramsey  fringes,  in  which  two  rr/2  microwave  pulses  are  applied  sequentially  to 
the  qubit.  Figure  6(a)  plots  the  SQUID  switching  probability  versus  the  time  interval  between  the  two 
pulses  for  a  microwave  frequency  fm  that  is  slightly  off-reso nance.  A  fit  to  these  oscillations  yields  a 
decoherence  time  T2*  of  7.1  ±  0.3  ns.  The  Ramsey  fringe  frequency  is  expected  to  be  proportional  to  )fm  - 
fq|,  where  fq  is  the  resonant  frequency  of  the  qubit.  This  behavior  is  clearly  demonstrated  in  Fig.  6(b),  which 
yields  fq  =  10.22GHz. 
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Qubit  flux  (0>o)  Microwave  amplitude  (a.u.) 

Fig.  7  (a)  Spectroscopic  peaks  at  varying  microwave  amplitude.  Excitation  frequency  was  10.2  GHz,  (b) 

Peak  linewidth  as  a  function  of  microwave  amplitude. 


Another  measure  of  decoherence  may  be  obtained  from  the  dependence  of  spectroscopic  linewidth 
on  microwave  amplitude.  In  the  limit  of  strong  driving,  this  dependence  is  approximately  linear  and 
extrapolates  at  zero  driving  amplitude  to  a  linewidth  corresponding  to  the  inhomogeneous  dephasing  time 
I2\  Figure  7  shows  linewidth  measurements  at  a  qubit  resonant  frequency  of  10,2  GHz,  and  the 
extrapolation  of  T2,=  l  1.1  ±  1  ns  from  this  data. 
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We  performed  echo-corrected  Ramsey  fringe  measurements  in  which  a  n  microwave  pulse  was 
inserted  between  the  two  k/2  pulses.  Data  for  a  range  of  n  pulse  positions  is  shown  in  Fig,  8(a).  The  echo 
pulse  corrects  for  inhomogeneous  dephasing,  so  that  the  peak  of  the  echo  envelope  decays  with  the 
homogeneous  dephasing  time  T2.  In  Fig.  8(b)  we  obtain  a  value  T2  =  30  ±  6  ns.  The  overall  dephasing  rate 
is  expected  to  be  the  sum  of  the  homogeneous  and  inhomogeneous  contributions,  that  is,  (Tj*)'1  =  (T2)'1  + 
(T2T1.  We  find  [(T2)'1  +  OVy1]"1  =  8,1  ±  1.7  ns,  which  is  consistent  with  the  measured  value  T2*  =  7.1  ± 
0.3  ns. 


IV.  Effect  of  Hot  Quasiparticles  in  the  Readout  SQUIB  on  the  Relaxation  Time  of  a 
Flux  Qubit 

The  flux  states  are  read  out  with  a  de  SQUID  inductively  coupled  to  the  qubits.  This  device  offers 
high  flux  sensitivity  readout  due  to  the  highly  nonlinear  dynamic  conductance  of  the  SQUID  when  operated 
at  a  bias  point  near  its  critical  current.  However  the  switching  of  die  SQUID  to  the  voltage  state  also 
results  in  the  breaking  of  Cooper  pairs,  which  compromises  the  low-noise,  non-dissipative  quality  that 
makes  superconducting  devices  promising  candidates  for  quantum  computing. 
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Fig.9  Qubit  relaxation  time  as  a  function  of  repetition  time  for  two  samples.  Fits  from  the  model  (solid 
red  and  blue  curves)  are  shown  overlaying  the  data  points.  The  open  fit  parameters  for  each  plot  are  given. 


The  number  density  of  these  nonequilibrium  quasiparticles,  and  hence  the  equivalent  Nyquist 
noise  spectrum,  is  roughly  proportional  to  the  SQUID  critical  current  and  the  length  of  time  the  SQUID 
remains  at  the  2A-gap.  The  quasiparticles  are  removed  through  pair-recombination,  a  process  that  takes 
place  on  a  time  scale  dependent  on  their  effective  temperature.  Recombination  lifetimes  in  Al  at  thermal 
equilibrium  with  the  bath  can  be  extremely  large  compared  with  nanosecond  timescales  required  of  a  qubit. 
For  pulsed  readout,  the  effect  of  the  long  recombination  tune  should  be  directly  observable  by  setting  a 
limit  on  the  repetition  rate  of  the  current  pulses  in  the  SQUID.  Below  this  limit,  the  rate  of  quasiparticle 
generation  exceeds  the  rare  of  quasiparticle  removal  in  the  SQUID,  and  qubit  decoherence  times  may  be 
drastically  reduced  by  the  resultant  increase  in  dissipative  noise  currents  in  the  SQUID,  which  appear  as 
flux  noise  in  the  qubits.  Indeed,  we  have  observed  [Fig.  9]  this  effect  in  relaxation  time  measurements  of 
two  separate  samples.  We  see  the  upper  limit  to  the  repetition  rate  in  the  relaxation  time  measurements  for 
both  samples  to  be  around  1  kHz.  A  viable  qubit  readout  device  will  require  a  much  greater  repetition  rate. 

To  understand  this  effect  and  to  explore  the  possibility  of  mitigating  the  limitations  it  imposes,  we 
have  developed  a  model  for  the  quasiparticle  generation  and  recombination  in  our  SQUID.  We  take  into 
account  Two  effects  of  the  high-energy  recombination  phonons  in  the  aluminum  film  of  the  SQUID:  the 
pair-breaking  process  and  phonon-trapping.  We  assume  the  quasiparticle  recombination  rate  per  unit 
density,  R ,  and  the  pair-breaking  lifetime  of  the  2A-phonons,  iB,  are  constants  and  use  typical  theoretical 
values.  We  also  assume  that  the  gap  anisotropy  smearing,  g,  is  constant  at  O.OSAaj,  The  phonon-trapping 
factor,  Fm,  is  the  average  number  of  pairs  a  given  phonon  will  break  before  escaping  to  the  substrate.  We 
leave  F^  as  a  free  parameter.  The  measured  background  relaxation  time  of  the  qubits  is  To,  which  is 
independent  of  the  repetition  rate,  Trtp‘\ 
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Fig.  10  Microwave  and  SQUID  current  pulse  sequence,  showing  the  resultant  nonequilibrium 
quasiparticle  density  in  the  read-out  SQUID.  The  qubit  samples  the  noisy  “hot”  quasiparticle 
environment  just  before  each  subsequent  microwave  pulse. 


For  every  point  in  a  train  of  pulse  sequences  [Fig.  1 0]  in  the  SQUID  and  qubit  system,  we  solve 
the  coupled  Roth  warf-Tay  lor  equations,  which  describe  the  quasiparticle  and  2A-phonon  densities  in  the 
SQUID.  The  effect  of  the  phonons  is  largely  to  retard  the  effective  recombination  rate  of  the  quasiparticles. 

The  quasiparticle  density  solution  has  three  regimes  of  interest:  the  driving  regime,  when  the 
SQUID  is  switched  and  generation  rate  exceeds  the  recombination  rate;  the  saturation  regime,  when  the 
recombination  rate  balances  the  generation  rate;  and  the  decaying  regime,  when  the  SQUID  is  restored  to 
the  zero-voltage  state  and  the  quasiparticle  densities  decrease  hyperbolically  with  time. 

The  qubit  is  turned  on  with  a  microwave  pulse  after  the  SQU  ID  has  been  in  the  decaying  regime 
for  a  time  -trep.  The  quasiparticle  number  density  is  calculated  at  this  time,  and  converted  to  an  effective 
temperature,  T*,  by  the  usual  BCS  relation.  We  then  calculate  the  equivalent  quasiparticie  noise  current 
spectrum  in  the  SQUID  ai  temperature  T*  and  subsequently  the  flux  noise  it  induces  on  the  qubits.  From 
the  standard  spin-boson  treatment,  we  arrive  at  a  qubit  relaxation  time  from  this  quasiparticie  number 
density.  By  repeating  the  calculations  for  many  repetition  times  in  the  range  of  25-5000  ns,  we  obtain  the 
solid  curves  shown  in  Fig.  9. 


By  solving  the  coupled  2A-phonon  and  quasiparticie  rate  equations  and  using  the  measured 
parameters  and  one  fit  parameter  (F^)  we  obtain  a  reliable  fit  for  the  relaxation  time  dependence  on  the 
current  pulse  repetition  rate  for  the  two  samples  at  two  different  frequencies.  It  is  particularly  gratifying 
that  the  best  fit  for  two  different  samples  is  obtained  for  the  same  value  of  F^.  From  this  analysis,  we 
conclude  that  efficient  phonon-substrate  coupling — particularly  around  the  SQUID  junctions — will  greatly 
increase  the  effective  quasiparricle  recombination  rates.  Additionally,  since  the  quasiparticles  are  generated 
in  the  junctions  the  presence  of  normal-metal  traps,  each  with  a  high-quality  NS  interface,  near  the  SQUID 
junctions  should  drastically  reduce  the  diffusion  of  large  quasiparticie  densities  around  the  SQUID  loop 
and  thus  minimize  the  resultant  flux  noise  in  the  qubits.  Alternatively,  one  should  adopt  an  inductive 
readout  scheme  in  which  the  SQUID  remains  in  the  zero  voltage  regime. 


V.  Entangling  Flux  Qubits  with  a  Bipolar  Dynamic  Inductance 

To  implement  a  quantum  algorithm,  one  must  be  able  to  entangle  multiple  qubits,  so  that  an 
interaction  term  is  required  in  the  Hamiltonian  describing  any  two  qubit  system  [6],  For  two 
superconducting  flux  qubits,  the  natural  interaction  is  between  their  magnetic  fluxes.  Thus,  for  two  flux 
qubits  arranged  so  that  a  flux  change  in  one  qubit  can  alter  the  flux  of  the  second  qubit,  the  coupled-qubit 
Hamiltonian  takes  the  form 


H  =  H]  +  H2  -  (K/2)ctz(1!  <jz(2)  . 


(1) 
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Here,  H]  and  H2  are  the  Hamiltonians  of  the  separated  qubits  1  and  2,  and  az(1)  and  az(2)  are  the  Pauli  spin- 
1/2  operators  for  qubits  1  and  2.  The  parameter  K  characterizes  the  coupling  strength;  for  K  <  0,  the 
minimum  energy  configuration  for  this  system  is  with  the  qubits  antiparallel.  For  two  qubits  coupled 
through  a  mutual  inductance  Mqq,  the  interaction  strength  K0  is  given  by 

K0  =  -2Mqq|Iq(1)||Iq(2)|,  (2) 

and  will  always  be  present  in  any  practical  layout  [7];  Iq(1)  and  lq,2i  are  the  supercurrents  in  qubits  1  and  2. 

Pulse  sequences  for  generating  entanglement  have  been  derived  for  several  superconducting  qubits 
with  fixed  interaction  energies  [8],  [9].  However,  entangling  operations  can  be  much  more  efficient  if  the 
interaction  can  be  varied  and  ideally  turned  off  during  parts  of  the  manipulation.  Although  switchable 
couplings  have  been  proposed  for  flux  qubits  [10],  [11],  these  approaches  require  the  rapid  switching  of 
fluxes  on  the  order  of  a  flux  quantum  O0  which  is  difficult  to  implement.  We  have  devised  a  coupling 
scheme  for  flux  qubits  whereby  the  interaction  is  adjusted  by  changing  a  relatively  small  current  which  can 
be  varied  rapidly  with  existing  pulse  technology.  For  suitable  device  parameters,  the  sign  of  the  coupling 
can  also  be  changed,  thus  making  it  possible  to  null  out  the  direct  interaction  between  the  flux  qubits. 
Furthermore,  the  same  device  can  be  used  to  vary  the  coupling  and  to  read  out  the  flux  states  of  the  qubits. 

A  variable  coupling  scheme  for  charge-based  superconducting  qubits  has  been  suggested  recently  [12],  but, 
to  our  knowledge,  no  feasible  variable  coupling  has  previously  been  proposed  to  exploit  the  advantages  of 
flux  qubits. 

The  coupling  mechanism  for  our  scheme  [Fig.  1 1(a)]  is  the  circulating  current  J  in  a  dc 
Superconducting  QUantum  Interference  Device  (SQUID),  in  the  zero  voltage  state,  which  is  coupled  to 
each  of  the  two  identical  qubits  through  an  identical  mutual  inductance  Mqs.  A  variation  in  the  flux  applied 
to  the  SQUID,  <3>s,  caused  by  one  of  the  qubits  switching  between  its  two  flux  states,  changes  J  as 
determined  by  LD_1  =  3J/3d>s,  the  inverse  dynamic  inductance.  This  change  in  J  alters  the  flux  coupled  from 
the  SQUID  to  the  second  qubit,  leading  to  an  interaction  energy  between  the  two  qubits,  Ks,  which  is 
proportional  to  LD"\  The  coupling  strength  Ks  takes  the  form 

Ks  =  -2Mqs2  |Iq(1)||Iq(2)|  Re(ar/d<Ds) .  (3) 
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Fig.  11.  (a)  Layout  of  SQUID-based  coupling  scheme.  Element  Y  represents  the 

admittance  of  the  SQUID  bias  and  measurement  circuitry,  (b)  Response  of  SQUID 
circulating  current  J  to  applied  flux  Os  for  yic(0.45  O0)  =  0,  0.4,  0.6,  0.85  and  (3L  =  0.092. 

Inset  in  lower  right  shows  variation  of  J  with  Ib  near  Os  =  0.45O0  and  upper  left  inset 
shows  modulation  of  Ic  with  Os. 

For  reasonable  device  parameters,  the  sum  of  Ks  and  the  direct  mutual  interaction  between  the 
qubits,  Ko,  can  be  as  large  as  -0.3  GHz.  For  a  particular  range  of  values  of  the  screening  parameter  pL  = 
2LIo/d>0  and  for  fixed  magnetic  flux,  LD‘ 1  can  be  reduced  to  zero  and  even  made  negative  by  changing  the 
bias  current  Ib  [Fig.  1 1(b)];  here  Ib  <  Ic(Os),  the  critical  current  for  which  the  SQUID  switches  out  of  the 
zero  voltage  state  at  zero  temperature  in  the  absence  of  quantum  tunneling.  The  ability  to  change  the  sign  of 
Ld_1  makes  it  possible  to  null  out  Ko,  yielding  a  vanishing  net  interaction  so  that  the  qubits  are  truly  isolated 
[Fig.  12].  In  previous  single  qubit  experiments  [4],  including  our  own  measurements,  the  qubit  state  is 
determined  by  coupling  the  qubit  screening  flux  to  a  similar  dc  SQUID  to  which  rapid  pulses  of  Ib  are 
applied  to  measure  Ic(Os,  T).  Thus,  the  existing  measurement  technology  allows  for  Ib  to  be  varied  on  this 
same  timescale  of  a  few  nanoseconds,  so  that  a  single  dc  SQUID  can  be  used  to  measure  the  flux  state  of 
two  qubits  and  to  couple  them  together  controllably. 


Fig. .  12.  Variation  of  K  with  Ib  for  Os  =  0.45  O0  and  device  parameters  given  in  text. 

When  combined  with  Rabi  pulses  to  generate  single  qubit  rotations,  our  coupling  scheme  can 
perform  entangling  operations  on  die  two  flux  qubits.  Our  calculations  of  the  coupler  dynamics  show  that 
the  quantum  Controlled-NOT  (CNOT)  logic  gate  between  two  flux  qubits  can  be  achieved  with  the  pulse 
sequence  shown  in  Fig.  13.  After  the  qubits  are  manipulated  individually  with  tuned  microwave  pulses, 
Si(t)  and  s2(t),  the  coupling  strength  K  is  pulsed  on  for  about  10  ns,  and  a  final  set  of  microwave  pulses  is 
applied.  The  total  duration  for  this  operation  is  about  30  ns.  Our  analysis  of  the  dissipation  in  the  bias 
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circuitry  indicates  that  the  decoherence  due  to  the  coupler,  as  characterized  by  the  Leggett  dissipation 
parameter  a,  can  be  reduced  to  about  10“*  for  a  reasonable  choice  of  series  resistance  on  the  bias  current 
leads,  2.4  kfl  This  value  of  a  corresponds  to  a  decoherence  time  of  500  ns,  which  is  comparable  with 
times  currently  achieved  experimentally  at  the  degeneracy  point,  but  significantly  longer  than  those 
observed  away  from  the  degeneracy  point. 


0  5  10  15  20  25  50 


t  (ns) 

Fig .  13.  Pulse  sequence  for  producing  CNOT  gate.  Bias  current  pulsed  to 
turn  on  interaction  in  central  region,  while  pulses  of  microwave  flux  drive  the 
qubit  biases,  £[(t),  e:(t),  for  producing  single  qubit  rotations. 

Finally,  we  note  that  this  approach  to  entangling  flux  qubits  should  be  readily  scalable,  for  example, 
with  the  layout  of  Fig.  14.  With  the  application  of  current  pulses  of  appropriate  magnitudes,  each  SQUID 
can  be  used  to  entangle  and  readout  two  qubits. 


Fig.  14.  Chain  of  flux  qubits  with  dc  SQUIDs  arranged  to  provide  variable  nearest  neighbor 
coupling  along  with  readout  of  qubits. 


VI*  Concluding  Remarks 

During  the  course  of  the  project,  we  have  developed  a  fabrication  process  for  flux  qubits  and 
assembled  the  necessary  infrastructure  to  characterize  them.  We  have  made  detailed  studies  of  the 
decoherence  times  obtained  from  Ramsey  fringes,  spectroscopic  line  width,  and  spin  echoes  and  shown  that 
they  are  self-consistent.  We  have  made  careful  measurements  of  the  effects  of  nonequilibrium 
quasiparticles  generated  in  the  readout  SQUID,  and  developed  a  theory  that  is  in  excellent  agreement  with 
the  data.  From  these  results,  we  conclude  that  it  is  imperative  to  develop  alternative  readout  schemes  that 
do  not  require  the  SQUID  to  switch  to  the  voltage  state,  for  example,  by  measuring  the  change  in  the 
inductance  of  the  SQUID  with  enclosed  magnetic  flux.  We  have  completed  a  theory  for  a  novel  technique 
to  entangle  qubits  using  the  readout  SQUID  as  an  inductive  entangling  device.  We  hope  to  implement  this 
scheme  in  the  near  future. 
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V 


We  did  not  succeed  in  our  stated  goal  of  entangling  two  qubits  experimentally.  This  failure  was 
due  entirely  to  the  virtual  demise  of  the  electron-beam  writing  system  in  the  Microfabrication  Facility  at 
UC  Berkeley  over  the  first  half  of  2005.  The  system  has  degenerated  to  the  point  where  it  is  almost 
impossible  to  achieve  consistent  results  from  run  to  run.  Despite  the  great  difficulties  in  writing  devices 
with  the  appropriate  parameters,  after  some  8  months  of  endeavor,  we  succeeded  in  producing  a  chip  with 
two  qubits  that  is  currently  under  investigation. 

Three  more  papers  are  currently  being  written  for  submission  to  journals  by  the  end  of  2005:  (a) 
A  paper  on  the  theory  of  the  three-junction  flux  qubit  in  the  limit  where  the  loop  inductance  is  non- 
negligible;  (b)  a  paper  describing  in  detail  our  experiments  on  quantum  coherence  in  the  flux  qubit;  (c)  a 
paper  describing  our  experiments  on  the  decoherence  in  the  flux  qubit  induced  by  nonequilibrium 
quasiparticles,  along  with  the  theory  of  these  effects. 
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